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Abstract 

Although the methods of calculation of power corrections in QCD sum rules are well known, 
algebraic complexity rapidly grows with the increase of vacuum condensates' dimensions. Cur- 
rently, state-of-the-art calculations include dimension 7 and 8 condensates. I summarize and 
extend algorithms of such calculations. First, I present all the formulae necessary for application 
of the systematic classification of bilinear quark condensates proposed earlier, and extend this 
method to the case of gluon condensates. Then I apply these systematic procedures to expan- 
sions of bilinear and noncollinear quark and gluon condensates in local ones, and of noncollinear 
condensates in bilocal ones. The formulae obtained can be used for calculation of correlators in- 
volving nonlocal condensates, and for inventing consistent anzatzs for these condensates. Finally, 
I briefly summarize the methods of calculation of heavy- and light-quark currents' correlators. 
This paper is aimed both to present new results on gluon and nonlocal condensates and to be a 
self-contained handbook of formulae necessary for calculation of power corrections in QCD sum 
rules. 
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1 Introduction 



The Operator Product Expansion (OPE) is widely used for investigation of correlators in quan- 
tum field theory. In particular, it is the basis of the QCD sum rules method. Let n(p) be a 
correlator of two local currents in the momentum space. Then OPE reads 

n(p)=5:c.(p2)(o,), (1.1) 

i 

where (Oj) are vacuum averages of local operators of various dimensions d (called also vacuum 
condensates), and Ci{p'^) are coefficient functions. They are perturbative series in a^. In order to 
make sum rules more precise and to establish their applicability regions it is necessary to calculate 
both higher perturbative corrections (higher terms in ag expansions of coefficient functions) and 
higher non-perturbative (power) corrections. I shall not discuss physical applications of QCD 
sum rules here. 

Methods of calculation of higher loop diagrams (perturbative corrections) form an established 
field of research, and many articles and reviews are devoted to them. I shall not discuss this 
subject; instead I shall concentrate on the methods of calculation of higher-dimensional terms 
in the OPE ( [L.l| ) in the leading order in Ug- 

There already exists an excellent review on technical problems of QCD sum rules calcula- 
tions 1^. The present article can be considered as an addendum to it. In all cases when a 
question is discussed in 0], I shall refer to this review and the literature cited in it. I should like 
to apologize to the authors of these articles for not citing them directly. 

I summarize and extend the methods of calculation of higher power corrections to correlators 
of heavy- |3| and the light -quark currents [4-9] . Full technical details (partly omitted from 
the original papers) will be included for the reader's convenience. In many cases, even if an idea 
of a calculation technique have appeared in the literature, it is very tedious and error-prone to 
rederive all the necessary formulae for each specific application. A set of reliable formulae can 
provide an invaluable help for many calculations of power corrections to QCD sum rules. 

I shall formulate the algorithms of calculation of higher power corrections in a form suitable 
for a Computer Algebra implementation. Unfortunately, there is no complete package for such 
calculations in conventional Computer Algebra systems (REDUCE, Mathematica, . . . ) . Such 
a package would be very useful. I have written several REDUCE procedures and Modula- 
2 programs (producing a REDUCE readable output) which implement some of the discussed 
algorithms. They are useful as a set of tools though they can't be a substitute for a complete 
package. 

The plan of the article is following. In Section |2| I briefly introduce the fixed-point (Fock- 
Schwinger) gauge on which the modern methods of calculation of power corrections are based; 
for more details and proofs see [|^ . Section ^ is devoted to the systematic classification of vacuum 
condensates. There exist many relations between condensates, therefore one has to choose a basis 
of independent condensates and formulate an algorithm of reducing an arbitrary condensate to 
this basis. A systematic classification of bilinear quark condensates has been proposed in P]. 



A carefully checked set of formulae implementing this method is presented in Section 3.1; this 



makes its application to any particular case very easy. In Section |3.2| , I extend this method to 



the case of gluon condensates. Section 3.3 contains some useful techniques of averaging two- 
and three-point correlators. 

Any calculation of power corrections to a correlator naturally involves nonlocal condensates. 
In the standard approach, they are expanded in series in local ones. Section ^ is devoted to 
a detailed study of some simple types of nonlocal condensates using the systematic methods 
introduced in Section |^. Tree-level calculations of two-point correlators usually involve gauge- 



invariant bilocal condensates. Their expansions in local ones are considered in Section 4.1. 



2 



Three-point correlators, as well as correlators involving bilocal currents (which are used in 
sum rules for hadron wave functions) require more complicated trilocal objects — noncollinear 
condensates (Section [4.2D . In some applications, such as sum rules for wave function moments, it 
is convenient to expand noncollinear condensates in series in one of the spacings. Coefficients in 



these series are bilocal condensates. Such expansions are constructed in Section 4.3. Expansions 
obtained in this Section can be used for calculation of correlators involving nonlocal condensates, 
and for inventing consistent anzatzs for these condensates. 

Finally, I briefly summarize methods of calculation of heavy quark currents' correlators 
(including heavy quark condensates) in Section ^, and of light quark currents' correlators — in 
Section |6[ 

2 Fixed— point gauge 

The QCD vacuum has a non-trivial structure. In order to calculate a correlator in QCD, we 
should first calculate it in a background (vacuum) gluon and quark field. Then we should 
average the expression for a correlator via these fields over the vacuum. After that we obtain 
the expression for a correlator via vacuum condensates ( |1.1| ). 

Correlators of colourless currents are gauge invariant. Therefore we can use any gauge for 
the background gluon field. It is convenient to use the fixed-point (Fock-Schwinger) gauge 

x^A^ix) = 0. (2.1) 

In this gauge the Taylor expansions for A^(a;) and q{x) can be written in a gauge-covariant 
form im 

A1{x) = ^x.Gl^ifi) + J-^x„x,D«G^^(0) 

+ -^^xpXc.x,DpD^Gl^{Q) + 0{x^), (2.2) 
q{x) = g(0) + Xe,Daq{Q) + ^x^x„D^D„g(0) + 0{x^) 

(and hence q{x) = g(0) + q{Q)DaXa + •••)■ Here D^q = {d^ — iA^)q is a covariant derivative 
in the fundamental representation, A^ = gAli^t"-; D^Gp^- = {d^i5"-^ + A'^^)G^p^ is a covariant 
derivative in the adjoint representation, A'^^ = gf°''^^A'i^. Hence only gauge covariant quantities 
are used at all stages of calculations. 

In this gauge the theory is not translational invariant. A translation should be accompanied 
by a gauge transformation to another fixed point gauge. Correlators become gauge invariant 
and hence translationally invariant after vacuum averaging. We can choose the gauge fixed point 
(origin in ( |2.1[) ) at any vertex of the correlator to simplify calculations. In the momentum space, 
background gluon lines have incoming momenta fcj and contribute A'^{ki) ^^^^'4 where ^^(fcj) is 
the Fourier transform of (|2.2D and is the series in derivatives of 5{ki). The vertex chosen as the 
gauge fixed point provides a common sink for all vacuum momenta fcj. 

The quark propagator (Fig. |l]) can be written as a sum over the number of background gluon 
lines: 

00 " _j_ 

^{p) = ^k{p), So{p) = ^— — 2 ' ^k{p) = -So{p)ASk-i{p), 

k=0 P ^ 

Afj, = — (/ + 2)/! ' ' ' ^"^'^"m^"*^"! ■ ■ - ^ciidv (2-3) 



3 



iS{p) 



P 



P 



iSo{p) 



P 



p-ki p - ki - k2 



+ 



+ ■ 



P 



ki 



p-ki 



+ 



Figure 1: Quark propagator in background gluon field 



It is convenient to use recurrent relations for Sk{p)- 



1=0 



oo < 

1 + 2' 



I 



:daDaAkl-i. 



(2.4) 



Here = d/dp^ acts on all propagators to the right, and — only on the nearest G 
9Ga,^{0)t°' ■ For example, retaining gluon operators with d < 4 (Fig. |l|) we have 



i 1 

S = Sq + -G^iySod^'ji^So + -DaG^iySodad^^^So 

i 1 

— -^DpDaGfj^ySQdpdad^'ypSQ — -Gpc^G^uSodp^aSod^^,ySo. 



(2.5) 



I have written a Modula-2 program for simplification of expressions constructed from 5*0, 
d^, and 7^. It applies to all 5*0 to the right using 5^5*0 = —Sq^^Sq, collects similar terms, 
and produces a REDUCE readable output, which can be used, for example, to expand these 
expressions in basis 7 matrices. 

These formulae can be used for both massive and massless quarks. For a light quark we can 
expand 5o in m. For a massless quark the equation (p^) gives 



S{p) = ^P- ^P,^Gf,^j^j5 



(2.6) 



+ 



1 



p J - p^J^p + p\Dxp^Gpy-fu) - 2ipxDxp^G^ulvl5 



2 {ph,G,xGx.P. - P,G,xGx.P.p) - 2^pxDx {p^J - p,J,p 



+ (A{pxDxf-p^D^) p^G 



P^iGpx.Gxvlvlb 



Here = gJ^t"-, = D,yG^j^ = gJ2q' (fltif^Q'- Note that the last term in ( |2i^ ) is missing in [| 
In deriving this formula, we have used the relations D'^G^y = —D^J^ + DyJ^ + 2i[Gp_x,Gxu 
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DyD\G pi^y — D\J^ — i[Gyx,G G ^xG + GyxGxfi — — 2^ ^uG puG , Gp^Gpy — Gp^Gpy. In 
practical calculations it is often more convenient to substitute symbolic expressions for propa- 
gators like ( |2.5D to the diagrams and to simplify the complete answers rather than to use ( |2.6| ). 

In order to consider the gluon propagation in the background gluon field we should substitute 
^d^ + a^ into the QCD lagrangian, where A'^ is the background field and is the quantum gluon 
field. The lagrangian becomes 

L = -\g%G% - \{Dpal){Dpal) + \{Dpal){D,al) + ^<G^ta^ + • • • (2.7) 

Here G^^, and Dp include only the background field yl^, and = gf'^^'^G'p^,. We choose 
the fixed-point gauge for the background field This lagrangian is still gauge invariant 

with respect to a°. We should add a gauge fixing term and a corresponding ghost term. It 
is convenient to use a generalization of the Feynman gauge fixing term — ^(D^aJ^)^. Then the 
free gluon propagator is D^^^{p) = Spu/p^ , and the vertices of the gluon's interaction with the 
background field (Fig. |) are iAf(jpi + p2)xSp,y - ^^^^(pi - P2)u - 'iAf{p2 - pi)p - iG''^^, and 

Aac Acbs: Aac Acb j_ Aac Acb 

"fJ-f -^a -^v ■ 




Figure 2: Vertices of gluon's interaction with background gluon field 




Figure 3: Gluon propagator in background gluon field 
Retaining terms with d < 4, we have the gluon propag ator (Fig. D 

Dap = -^^al3 + ■^'^Gap + ^ (^iPfiJfiSaP + ^ipxDxGa/S 



1 

p8 



2pxDxPi.Ji.Sap - 2 (^ipxDxf - p'D^) G. 



aP 



(2. 



+ 2 [P^G^iyGf.u + "iPfiG^xGxuPu) Sai3 + Ap^GaxGxp ■ 
Here the matrix notations are used for colour indices. 



3 Vacuum condensates' classification 

Vacuum condensates can be divided into classes according to the number of quark fields in them. 
Those without quark fields are gluon condensates. Those with two quark fields are bilinear quark 
condensates; they have d > 3. Four-quark condensates have d > 6, and so on. The unit operator 
is the gluon operator with d = 0, according to this classification. 

As is clear from the previous Section, bilinear quark condensates of the form 
{q{D . . . DG){D . . . DG) . . . 7 . . . ^q) appear in calculations of bilinear quark currents' correla- 
tors at the tree level. Analogously, gluon condensates of the form (Tr(L' . . . DG){D . . . DG) . . .) 
appear in the one-loop approximation (there are at least two (D . . . DG) groups because other- 
wise the colour trace vanishes). 



In Section 3.1 I discuss the systematic classification of bilinear quark condensates follow- 
ing 1^. All formulae necessary for the practical use of this method are presented. In Section 
I apply similar methods to the gluon condensates. Techniques of vacuum averaging of expressions 



for correlators via quark and gluon fields are discussed in Section 3.3 



3.1 Bilinear quark condensates 

Using the formulae G^^, = gG^^^f" = i[D^,D^], (D^A) = {DfA^)t'' = [D^,A] (where A'' is 
in the adjoint representation of the colour group, and A = A^-f^), we can easily reduce any 
bilinear quark condensate of dimension d = m + 3 to a combination of terms of the form 
{qD^^ ■ ■ ■ D^,^Tfj_^,„f,^q), where r^i...;^^ is constructed from 7^ and 5^,,. Choosing the terms 
with the largest number of 7 matrices, we permute them in such a way that their indices are in 
the same order as in Z)^^ . . . D^^ . Arising additional terms have fewer 7 matrices. We repeat 
this procedure going to terms with fewer 7 matrices until we reach terms with one 7 matrix or 
without them at all. As a result, a bilinear quark condensate reduces to a linear combination 
of terms of the form {qOiq), where Oj are constructed from D and D^. Due to the equations of 
motion, those terms in which D is adjacent to (7 or q reduce to lower-dimensional condensates 
multiplied by m. 

Having written down all the condensates Bf = (qOf^q) where Of* are constructed from D 
and D^ and no D is adjacent to q and q, we obtain a certain set of d-dimensional bilinear quark 
condensates. We have just demonstrated that any d-dimensional condensate can be systemati- 
cally expressed via Bf, mBf-\... This means that they form a basis of quark condensates with 
dimensions < d. But this basis is extremely inconvenient for use because its condensates contain 
a maximum number of derivatives acting on q. So we should better choose a basis of the most 
convenient d-dimensional condensates Qj. We can express Qj via Bf, mBf~\. . . Solving this 
linear system, we obtain the expressions for Bf via Qj, mQ'^~^,. . . Having these expressions, we 
can easily reduce any d-dimensional bilinear quark condensate to the basis ones 

A general guideline for choosing good Q'^ is to have a minimum number of derivatives acting 
on the quark field, and hence a maximum dimension of gluon operators in Q^. It allows to expand 
propagators to a minimum dimension in the background gluon field. Among such operators one 
should first of all choose those containing because they are really four-quark ones and are 
more easily calculable. They may be suppressed in some vacuum models (of the instanton type) 
in which vacuum gluon fields may be strong but D^G^u ~ 0. 
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For dimensions d < 6 we have 



B' = {qq), B' = -(^qD'q), B^ = -ii^qD^DD^q) . (3.1) 
We choose the basis condensates 

Q^ = {m), Q' =i{qG^,a^,q), = IqJq) : (3.2) 



5 TO ncnollTT r\ dirt nA- Ckr\ nn^/^^^ 



where a^i, = 7[^7i/] , and square brackets mean antisymmetrization (Q^ is usually denoted mgQ 



These condensates are expressed via Bf as 



q5 = _2B^ + 2m25^ (3.3) 
q6 = _2s6 ^ 2mB^. 



Solving this linear system we obtain 

B"^ = Q^ 



S^ = -^Q^ + "iV, (3.4) 



B[ = {qD'D'q) , i?^ = {qD^D^D^q) , 



For d = 7 we have 



S| = {qD^D^D^D^q) , fij = [qD^D^D^q) . (3.5) 

We choose basis condensates 

<3i = {qG^^G^uq) , Ql = i (gG^i.G^i.75g) 

= {qGfixGx,y<7n^q) , = i {qDixJuCr^uq) ■ (3.6) 

Here the condensate containing G^i, = ^e^^p^Gpa and 75 = ■^EaP'ySla^fSl'y^s is understood as a 
short notation for the expression from which both e tensors are eliminated using £^^'''^ £aj3'y& = 
— 4!(5j^(5^5;^(5^j , and this expression is valid at any space dimension D. These condensates are 
expressed via Bf as 

Q\ = 2^1 - 2i?|, 

Ql = -2BI - 2BI + 2BI + 2BI - AmB^ + 6m^B^ - 2m^S^ 

Ql = 2Bl - 2Bl -Bl + 2mB^ - rri^B^, (3.7) 

Ql = -2Bl- 2Bl + 45| - 4.mB^ + Arri^B^. 

Solving this linear system we obtain 



Bl 




-\q\ 


-Ql- m'^Q' 


' + m^Q'^, 


Bl 




-\qI 


-Ql + \q\ 


— mQ^ — 1 


Bl 


= Ql- 


\qI- 




mQ^ — m 


Bl 


= Ql- 


Ql - mQ' - ^Q' 


+ m^Q^. 



(3.8) 
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Similarly, for d = 8 we have 

Bf = i (liiD^DD^D^ + D^D^DD^)q] 
Bl = i (q{D^DD^D^ + D^D^bD^)q) 

Bl = i (q{D^DD,D^D, + D,D^D,DD^)q) , (3.9) 



Bl = i (qD^DD^q) , BI = i (qD^D,DD,D^q 
Bl = i (qD^D,DD^D,q) , B^ = i (qD^^D^D^q) . 

Operators similar to those in Bf 3 but with a minus sign between two terms are C-odd, and 
their vacuum averages vanish. 

An interesting new phenomenon similar to the axial anomaly arises for bilinear quark con- 
densates of dimensions d> 8. At a space dimension D 7^ 4 we can construct the condensate 



A = i (qDaDfsD^DsDs-f[aJf3jj75l£]q) ■ (3.10) 



All coefficients in the expansion of this condensate in gluon condensates (see Section 5.1) contain 
traces vanishing in the 4-dimensional space. Momentum integrals in coefficients at d > 8 gluon 
condensates converge, and at d = 8 they contain divergencies 1/e. As a result, the quark 
condensate A is equal to a combination d = 8 gluon condensates. 

Choosing the anomalous condensate A as one of the basis ones, we can choose 6 proper quark 
d = 8 basis condensates as 

Qf=i{q[[G^X,Gx,h,D^]+-f,q), 

Q2 = - {qiiGt^x, Gai.], I^^l+Ti^Tsg) , 

Ql=i(q[DG^,,G^,]q), Ql = (qD^Jq), (3.11) 

Here again the expressions containing G^j^y and 7^75 = 7^7/377 are understood as a short 

notation for the expressions in which both e tensors are eliminated. 

All operators in ( ^.llD (as well as in ( |3.6D , (|3.2|) ) are purely C-even. C-conjugation permutes 
all and 7^ in the opposite order, and changes their signs (because is transformed to D 
Therefore every commutator gives a factor —1 to the C~parity, and every anticommutator gives 
+1; we should not use ordinary products if we want to obtain operators with a definite C-parity. 
According to these rules, G^y and all its derivatives D\G^y,... (including J^) are C-odd; 7^ 
and a^y are C-odd, while 1, 75, and 7^75 are C-even. 

The condensates ( p.lOD , ( 3.11| ) are expressed via ( p.9[ ) as 



A = -Bf - bI + bI + Bl + Bl- bI + Bl 

+ 2mBl + 2mBl - 2mSj - 2m5j + Zm^B^ - 4m^B^ + rrV'B^, 

4, 



Ql = Bf-Bl + 2BI - 2BI - 2mBl + 2mSj 



q8 = 2Bf - 2Bl - 2Bl + 2BI - 2S| + 25^ 

- 2mBl + 2mBl + 2mSj - 2m5j - 2m^B^ + 2m^B^, (3.12) 
Ql = -ABI + ABI + 4mBl - 4mBl, 
Ql = 2Bf + 2BI - 2BI - ABI - 2mBl, 
Ql = Bf + Bl- 2BI - 2mBl - 2mBl + AmB'J 
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Ql = 3B^ + ^2 - 2-B| - 2B: 



4 



AmBl - 2mBl + 4mSj - Am^B^ + Am^B^. 
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Solving this linear system we obtain 

bI = a + \q\-\qI + \qI-\qI + qI 

+ rnQl - ^Ql - 2mQl - m^Q^ - 2m^Q^ + 2m^Q^ 

bI = a-IqI-IqI-\qI-IqI + qI 

+ 3mQl - ^Ql - 2mQl + mQl - Sm'^Q^ - 2m^Q^ + 2m^Q^ 

bI = a-IqI-qI + qI 

+ 2mQl - ^Ql - 2mQl + mQl - Sm'^Q^ - 2m^Q^ + 2m^Q^ 
Bl = A + ig? - ^Ql + ^Ql + ^Ql + "^Ql - mQl - m'Q'> + m^'Q^ (3.13) 

Bl = \a- \q\ - \qI - \qI - \qI - \qI + \qI 

+ ^mQl - "^Ql - mQl + ^QI - 2m^Q^ - m^Q^ + m^Q^ 

Bl = \a- \q\ - \qI + \qI - \qI - \qI + \qI 

+ mQl - jQl - mQl + yQj - 2m^Q^ - m^Q^ + m^Q^ 
B'r=A-lQl + lQl 

+ mQl + - mQl - \m^Q^ - + ^5q3^ 

I have implemented the discussed algorithm as a Modula-2 program. It accepts an arbitrary 
bilinear quark condensate and expresses it via Bf producing a REDUCE readable output, which 
can be used to express the condensate via using the formulae (|3.4|) , (|3.8|) , (|3.13 ). 

The problem of four-quark condensates' classification is much more difficult. Even at 

d = 6 there exist infinitely many condensates ■ ■ ■ l^M^]q) (^Tmi • • • T^n]?), (^i^Tl^i • • • T^n]?) 

(qt°"jf_ii ■ ■ ■ T^n]?^- Those with n > 4 are anomalous and can be eliminated. 

3.2 Gluon condensates 

In this Section we apply similar methods to the classification of gluon condensates. For each 
dimension d = 2n we introduce a linear space. We shall call it the extended space. It is 
generated by the basis of formal sequences Ef constructed from in which each index is 
contained twice. Sequences obtained from each other by renaming indices, cyclic permutations 
(trace cyclicity), and reversing (C-parity) are considered equivalent. Physical gluon condensates 
can be systematically expressed via this basis, and form a subspace of the extended space. We 
can choose a basis in this subspace composed from the most convenient condensates Gj, and 
extend it to a basis in the whole space using a subset of Ef. After that we can expand any 
condensate in Ef and then reexpress it via Gj and the selected Ef. In fact Ef should not appear 
if we have indeed selected a complete basis in the physical subspace. 

General guidelines for choosing good basis condensates are similar to the quark case. First 
of all, we should include all independent condensates containing J^. The number of derivatives 
should be kept minimum. 

There is only one gluon condensate with d = 4, namely G^ = {Ti G ^^uG ^j^y) . The extended 
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space is 2-dimensional: Ef = D^D^, = D^D^D^D^; = 2Ef - 2E^. But all this is of no 
use here. 

At (i = 6 the extended space is 5-dimensional: 

Ef = D^D^D^, El = D^D,,D,D^D,, EI = D^D^D^D^, 

El = DxD^DxD.D^D,, Ef = D^D^D^D^D^D,. (3.14) 

There are 2 linearly independent physical condensates: 

G\ = i (Tr Gx^.G^,G,x) , = (Tr J^) . (3.15) 
They are expressed via Ef as 

Gl = Ef-^El + 3El- El, Gl = -2Ef + 8El - 2El - 4El (3.16) 
We can exclude £'4 . 5: 

1716 ^ 7716 I o rp6 ^ 7716 

El = -Gf - ^Gt - \eI + 3i?f - \eI (3.17) 

After that we can expand any d = 6 condensate in ( |3.14 ) and express it via (|t|) and £f 3; 
the last terms should not appear. I have verified this for all d = 6 condensates. 
At d = 8 the extended space is 17~dimensional: 

^8 = D'^D^D'^D'^, El = D^D'^D^D^D^Dy, 

El = D^D^D^D^D^, Ef = D^DxD^DxD^D^D,, 

El = D^DxD^DxD^D^, Ef = D^DxD^D.DxD^D,, 

El = D^DxD^D,DxD,D^, eI = D^DxD^D,D^D,Dx, 

El = D^D^D^D^D^D^, Ef^ = D^D^D^D^D^D^, (3.18) 

Ef^ = D^D,D^D,DpD„DpD„, Ef^ = D^D,D^DpD,D^DpD„, 

Ef-i = D^D^D^DpD^D^DpD,, Ef^ = DpD,DpDpD^D,D^Dp, 

Ef^ = DpD^DpDpD^D^DpD^, Ef^ = DpD^DpDpD^DpD^D^, 

Efj = DpD,DpD„DpD,DpD,. 

A set of independent d = 8 gluon condensates was found in We choose one condensate in a 
different way Q: 

Gi = {Tr GpiyGpuGajsGajj) ; ^2 = {Tr GpuGajjGpi/GalB) , 
Gs = {Tt GpaGauGiyjsGf^p) , G| = {Tl GpaGauGppGpiy) , 

Gl = t (Tr JpGp,J,) , Gl = i (Tr ,h[DxGp,, Gp,]) , (3.19) 
Gl = {TTjpD^Jp 

All condensates are written in an explicitly C-even form (even number of commutators) . They 
are expressed via Ef as 

G\ = -8El + AEl^ + AEl^, 

Gl = -8El,+AEl, + ^El,, 

Gl = Ef - 4El + AEf + 2El - AEf^ + Ef^, 

Gl = 2El - 4El + El + 2Ef^ - 4Ef^ + 2^184 + Ef^ , (3.20) 
Gl = Ef- 5EI + 2EI + AEl - 2EI + AEf - AEI - eI + Ef^ + AEf-^ - AEf^, 
Gl = -4EI + 4EI + 8EI - 8EI - 12EI + AEf^ + 8Ef-^ - 8Ef.;i + 8-^1^4, 
Gl = -2Ef + 8^;^ - 6EI - 8EI + 12^;^ - IGEf + AEI + 8Ef4. 
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We can exclude Ef^ ^j: 



El, 

IT'S 

^12 



-G\ + 2El 



E 



10! 



A ^ 4 



+ -El--El + -El + El 



-El + E'r + El 



-E'o 



E- 



13 

4 

^14 



7^1 ~ o'-^6 + o^7 



^i^f + \eI + i^f 



^8 



:El 



-E 



10) 



-El + 2ii;« 



9) 



\gUgu\gI-\ 



-El + -El 
1 



2eI + 2£;| + 2^;^ 



10' 



16 



-G^ + GI-GI- E^ + 2El - 2EI + AE^ - E, 



E 
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2 

lG? + iG^-G8 + 2Gl + 2G^ 



Si 



10! 



(3.21) 



1 

2"^ ' 4"^ " ' — ' — "° ' 2 
- 2£;f + 2,eI - 2EI + 4£;| - 2EI. 

After that we can expand any d = 8 condensate in ( 3.1^ ) and express it via ( |3.19| ) and Ef 
the last terms should not appear. 

I have written a Modula-2 program that expands any gluon condensate in Ef and produces 
a REDUCE readable output, which can be used for expressing this condensate via Gj and 
selected Ef. Using this program, I have verified that all (i = 8 gluon condensates indeed can be 
expanded in the basis ( 3.19| ). 



3.3 Vacuum averaging 

After obtaining an expression for a correlator via background fields we should average it over 
the vacuum. The first possible way is to write down the most general expressions for vacuum 
averages with free indices and to find unknown coefficients by solving linear systems. 
For the quark condensates with d < 6 we have 

/- \ - ^1 

\QaQp) — 22 P^^ 



{la, 



pa, 



(q^DaDfsqp) = 



Q ( + 3'^a/3 



2m^Q^5^p 



(3.22) 



J per 



{q^DaD/sD^qp) = 



pa 



One can easily obtain similar formulae with GafB and its derivatives by antisymmetrization over 
corresponding indices. 
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For the gluon condensates with d < 6 we have 



qA 



pa 



+ 



NCfD{D - 1)(L>2 _ 4) 
-{D + 2)GU{{Sap6,a6f,p -ifi^ i^)) -{p^ a)) - (q ^ /?)) 
{D - 4)G? + {D- 2)Gl\ {{{6^p6,a6f3p - ^ z^)) - {p ^ + (" ^ /?)) 
A{D - + 2{D - 2)Gt\ 6af3{6ppS,^ - {p ^ a)) 



(3.23) 



2Gf 



y ^pu^af^^pa/ - N2CpD{D -l)iD- 2] 



f 



abc 



We can substitute these relations to expressions for correlators and obtain scalar expressions. 
But these formulae become very complicated in higher dimensions. 

An alternative and simpler way is to average an expression for a correlator over p directions. 
This is done using the recurrent relation (n is even) 



P 



D + n-2 
or the explicit formula 



PpiPp2 • • • Ppn 



(p2)«/2 



• PMu + • • • + ^PlPr.Pp2 ■ ■ ■ Pl^n^l ) (3-24) 



D{D + 2)...{D + n-2] 



{^P1P2^P3P4 • • • ^Pn-iPn + 



(3.25) 



where the sum is taken over all (n — 1)!! methods of paring n indices. Here D is the space di- 
mension; the formulae can be used at any dimension including non-integer one (the dimensional 
regularization) . After such an averaging we obtain scalar vacuum condensates for which the 
methods of the previous Sections can be directly used. 

In calculations of three-point correlators we encounter the averages P^i • • • Ppm Qui ■ ■ ■ Qun over 
orientations of the pair p, q with a fixed relative orientation. These averages can be calculated 
using the decomposition q = {qp)p/p^ + q^, q\ = {p^q^ — {pqy)/p'^- First we average over 
orientations of q± orthogonal to p in {D — 1) -dimensional space using ( |3.25| ) with 5f^i, — > 5±pi, = 
^pu ~ PpPu/p^ ■ Then we average over orientations of p. It is clear that if there is only one vector 
q then the averages are given by the formula ( 3.25| ) in which one p^ is replaced by pq. Some less 
trivial averages are 



1 



PaPpqpQy 



{{D + l)p\^ -2{pqf)5o.p5, 



{D-l)D{D + 2) 

+ {D{pqf -p^q^){5ap8i3y + Sav^Pp 

pq 



pu 



PaP/sp-rqxqpQu 



{D-l)D{D + 2){D + 4) 
X '{{D + l)/g2 _ 2{pqf){6af3Sxp6^, + ■■■) 

+ {{D + 2){pqf - 3p\^){6ax6f3p6^, + • • • 



(3.26) 



PaPfSP-yPsqpQu 



{D-l)D{D + 2){D + 4) 
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{{D + 3)p'q' - i{pqy)6^u{SapS^S + 



+ {D{pqr - P q){5c,p5^^ikv + ■■■) 

Of course, these formulae can be also obtained by writing down the most general forms with 
unknown coefficients and solving linear systems. 



4 Nonlocal condensates 



As a sample application we consider expansion of nonlocal vacuum condensates in local ones. 
Nonlocal condensates naturally appear in calculations of correlators [10-14]; their expansion in 
spacings lead to the usual OPE (1.1). Expansions derived in this Section can be used either for 
easy obtaining of OPE of correlators involving nonlocal condensates, or for inventing consistent 
anzatzs for nonlocal condensates. 



4.1 Bilocal condensates 

The gauge-invariant bilocal quark condensate 12] has the structure 



{qME{^,x)q,{x)) 



im) 



fs{x 



pa 



E{0,x) = Pexp -i J A^{y)dyf, 



In the fixed-point gauge E{0,x) = 1, and 

{q{0)q{x)) {q (l + ^XaXf3DaDi3 + ^XaXfSX^XsDaDfsD^Ds + 



fs{x' 



{qq) 



{qq) 



tX + 



2\DB^ A\D{D + 2)B^ 



-X 



3Q? - iQl - 3Ql + Ql- 2mQ^ - 3m^Q5 ^ ^^AqS ^ 



4!L»(L> + 2)Q3 



fvix 



2 iD {q{0)xq{x)) 



(qq) 



iD Iqx ixaDa + ^XaXpX^DaDjjD^ + ^XaXpX^XsXeDaDjjD^DsDe H j q . 



{qq) x'^ 



B^ + 2mB^ 



m 



+ 



3!(Z? + 2)S3 

Bf + Bl + Bl + Bl + Bl + Bl + 2m{Bl + BJ, + S|) 4 



51{D + 2){D + 4)B^ 



x^ + 0{x^ 



m + 



3!(Z) + 2)g3 

bA - \Ql + \Ql - \Ql - 3Ql + hQl + 5m(3gl -QI- 3QI + Ql) 



5\{D + 2){D + 4)Q-^ 



(4.1) 



(4.2) 
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On the first step we use the averaging over x directions ( |3.25| ) and the definitions of Bf (|3.1| ), 
(|3^), and on the second step — the expressions for Bf via Qj (^74|), ( ^l8[ ), ( |3.13| ). 
The gauge-invariant bifocal gluon condensate |13, 14 1 has the structure 



pa 



D{D-l) 



) = D{D- 1) (TrG^,(0)Gp,(x)) /G^ 



dDi{x^) 



(4.3) 



+ 



where £^"''(0, x) is the P~exponent in the adjoint representation equai to S"'^ in the fixed-point 
gauge. We obtain the hnear system 



Dix^) + D,ix^) + ^D[ix^) 
D{x^) + Di{x'^) + x'^D[{x'^) = 



B{x^) 



where 



B{x^) 



{TvG^,{Q)G^,{x))=G^ + 



^^^{TrG,mGuxix)) 



Aix^ 
2D 

= G^ 



+ 



A2X' 



4\D{D + 2) 
Bix^ 



+ 



B2X^ 



(4.4) 



(4.5) 



2(L> + 2) 4!(L> + 2)(L> + 4) 



+ 0{x^ 



Here 



Tr Gfj^pD'^G I , 



Si = + (TrG^A(I?|.i?. + D,D^)G,x) , 
^2 = (TrG^,(Z)2z)2 + D^D'D^ + Z),i?^i?„Z)^)G^,/ , 

^2 = ^2 + ( Tr G^x{D,,D,D^ + D^D^D^ + D^D^D, + D.D^D^ 



Using the method described in Section 3^, we obtain 

A, = -4G? - 2Gl 
Bi = - AGl 



(4.6) 



A-, = 3G?_2 + 24G^_4 - 36G^ + lOG^ - 

5G?_2 + 50G^_4 - 84G^ + 24G^ - 18G?, 
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B2 



where Gf_2 = G? - G^, Gi_4 
D{x^ 



Gi 



G?. The solution of 



is 



1 

G4 



(2^1 - Bi)x^ (3^2 - ^2)^^ 



2{D - 2) 



+ 



+ 0(x^ 



Gf 



(L> - 2)G4 

= 1+ ' " 



2 I 2^1-2 + 11^*3-4 



2-4!(D2 -4) 

12G| + 3G| 



G^D 



(L»2 _ 4)^4 

Z? + 2 D 



4! 



+ 0(x6), 



(4.7) 



A. 



X 



Ai\ x' / B2 

2 \D + 4 D J 2-4!(L» + 2) 



+ 0{x^ 
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+ 



{D - 6)G?^2 + (13-D - 48)G^_4 - 24(Z) - 3)G^ + {ID - 20)G^ 

,4 



6{D - 2)G 



X 



4!(D + 4)(L>2 -4)G4 



0{x'° 



In the abelian case, only Di{x^) survives ||T^; this is the reason why D{x'^) ( |4.7| ) contains only 
purely nonabelian condensates G\, G\_2-, G3_4, G|g. 



4.2 Noncollinear condensates 

Three-point correlators, sum rules for hadron wave functions, and some other kinds of cal- 
culations naturally involve trilocal objects which we shall call noncollinear condensates. The 
noncollinear quark condensate has the structure 



(gJy)E(y,0)£;(0,x)gp(x)) 



(4.8) 



fsix,y) + 



1 



J pa 



In the fixed-point gauge E{y,0) = E{0,x) = 1. The difference 
fs{x, y) - fs{{x -y) ) = 



A {q{2D^D^ - D^D^D^ - D^D,D^D,)q) 



+ 0{x\y^) 



12D{D - 1) {qq) 



(4.9) 



A (3Ql + 2Ql - AmQ^ 
2AD{D - l)g3 



3 „,3\ 



0{x\y 



(where A = x^y^ — (xy)^) vanishes when x and y are collinear, as expected. The coefficient 
/t^x, y) multiplies the structure [x, y] vanishing in the collinear limit, therefore it need not vanish 
itself: 



fT{x,y) 



DiD-l){q{y)[y,x]q{x)) 



1 



(4.10) 



^ {6Ql - 3Ql - 6Ql + 2Ql - 3mQ^ - 3m^Q^) {x - yf - (3Qg + 2Ql) xy ^ 



6{D + 2)Q5 



For the vector matrix elements we find 



where 



i {q{y)xq{x)) = -i {q{y)yq{x)) 
Q^{x^ - xy) 



D 



-fviix-yY) 



A 



12D{D - 1 



-gix,y) 



(4.11) 



gix,y)=Q^ + 



1 



20A{x'^ - y^) - 30Q^(x^ - xy) + l(^Ql{x^ - y 



40{D + 2) 

- Ql{13x^ -y^- 9xy) - 2QI{2x^ + y^ - 6xy) - 2QI{7x^ + 6y^ - Uxy) 
+ 20Ql{y'^ - xy) + 60mQ^(x^ - xy) - 10mQl{x^ - y^) + AOmQl{x^ - xy) 

- 20m2Q6(5x2 + y^ - 6xy)l + 0(x^ y^). (4.12) 
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Solving the system of two linear equations, we finally obtain 

y'^9{x,y) +xyg{y,x) 



fv{x,y) = fv{{x-y) ) 



12{D - 1)Q3 

The noncollinear gluon condensate has the structure 



(4.13) 



GUy)E''\y,0)E'-iO,x)G%ix 



be/ 



K„ 



ni/yyj---^ \tii"j-^ pav^j / Df^D — 1) ^^^^''^'^ 

A 



{x,y), 



'-^nfpcr — G {^Kni/^p(j(x,y^ K^^ii^p(j(x yy) 

+ C2{{x^,Xij5ua - ^ i^)) -{p ^ cr)) 

+ c^{{y^iypKa - (m ^ T^)) - (p ^ o-)) 
+ Ci{{x^,yp5ya- - (/w ^ i^)) -{p^ o-)) 

+ C5{{y^Xp6ua - (/W ^ ^)) - (p ^ cr)) 

+ ceix^y^ - (p ^ i^)){xpya - (p ^ o-)) 
with C3{x,y) = C2{y,x). We find 



(4.14) 



A 



(-3Gt2 + 2Gi), 



Txn^Xi/Xfj^Xi/ = {T\^^\^y^y 
A 



(4.15) 



12(D + 2) 



(4^2 + xy) (-Gf _2 + Gi) - {x^ - 2xy)Gl_^ - 4xy Cl 
A 



12{D + 2) 



12{D + 2)G5' - (2x" + 2y^ + .Ty)Gf_2 



(llx^ + lly^ - 20xy)Gi_4 + 12(x - y)^G§ - (3x' + 3y^ - 10xy)G^ 



T\^^\iiy^Xi/ 



A 



12(Z) + 2)G5 + (x2 + y2_7^y)G'f_2 



i2(i:> + 2) 

+ (13x2 + 13y2 - 22a;i/)G^_4 - 8(2^2 + 2y'^ - 3a;y)G| + 2(2a;2 + 2y^ - xy)Gl 

A2 



Tfxii^paX fxX pyiiycr 



(-5Gt2+4Gi_4 + 12G«). 



12(L> + l)(L> + 2) 
Solving the system of 6 linear equations, we finally obtain 

A [-D{?,D - 7)Gf_2 + 4(L> + 3)G|_4 + 2{D -2){D- 2,)Gf\ 



Cl 



C2 



12(L> + 1)(£» + 2){D - 2){D - 3) 
^ {A{D + l){D-?,)xyGl 



12{D + 1){D + 2){D - 2){D - 3) 
+ {D-Z) [D{2y^ - xy) - Ay^ - xy] G^ 

(L»2 _ 3)(y2 ^ 2xy) + 2L»(3y2 - 2xy)] Gt4}, 
C3 = C4(x ^ y), 

^2A{D + 2)G? - 3(x - y)'Gt2 



(4.16) 



C4 - C5 



12(1)2 _ 4) 

6(4^2 + 4y2 - 7xy)G|_4 + (Tx^ + Ty^ - 12xy)(4G^ - Gl) 
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C4 + C5 = ; TT \-, TT r( - 4(D + 1)(Z) - 3)(x^ + y^)G\ 



12(L> + 1)(D + 2)(L» - 2)(Z) - 3) 
+ (Z) - 3) + - 4xy) + + + Sxy 



p8 



- yy - 2D{x^ + y - 6xy) - 3(x + yY Gl^2 



+ 2 



(L>2 - 3)(x2 + + xy) - 2D{x'^ + y"^ - 3xy) 



^3-4 f' 



_ (D - 4)Gf_, + - ifGU + iD- 2)P - 3)Gi 
2{D + l){D + 2){D-2){D-3) 

In the abelian case K^u,p<T{x,y) = K^u,p(t{x — y) has only the Di{{x — y)^) term (^^), therefore 
Tfj,u,pa contains only purely nonabelian condensates. 

All calculations in this Section were performed using REDUCE (with some REDUCE code 
generated by the Modula-2 programs described in Section |3|). Some of them require a lot of 
memory: I could not complete the gluonic calculations on a PC with 8 megabytes. 

4.3 Expansion in bilocal condensates 

In some applications, such as sum rules for wave function moments, it is convenient to expand 
noncollinear condensates in series in x keeping y finite. Coefficients in these series are bilocal 
condensates (Section ^ ). 

Expanding a noncollinear quark condensate {q{y)E(y,0)TE{0, x)q{x)) in x, we obtain terms 
like {q(y)E(y,0)Oq{0)) where O involves derivatives and gluon fields at a; = 0. Therefore we 



encounter the problem of classification of such bilocal condensates. Similarly to Section 3.1, we 
can easily reduce any condensate of this kind to a linear combination of terms with O constructed 
from D, {yD), and D^, maybe multiplied by y from the left. Terms with D at the right reduce 
to lower-dimensional condensates multiplied by m. Terms with (yD) at the left reduce to 
derivatives of lower-dimensional condensates: expanding the equality {q{y)E{y,0){Oq)x=ay) = 
(g((l — a)y)£^((l — a)y, O)(0q)x=o) in a we obtain 

{q{y)Eiy, 0){yD)Oq{0)) = " {q{y)E{y, O)Og(O)) . (4.17) 

There are 2 basis dimension 3 bilocal condensates: 

Bf{y^) = {q{y)E{y, O)g(O)) , Bl{y^) = -i {q{y)E{y, {))yq{Q)) . (4.18) 



They can be simply expressed via 2 bilocal functions considered in Section 4T 

Bl{y'') = Q'fsiy'), Bliy^) = Q^fviy')^. (4.19) 

Expansions of these functions at small y^ are given by equation ([4.2|). 
There are 4 basis dimension 5 bilocal condensates: 

B!{y^) = -{q{y)Eiy,0)D^q{0)] 
Bl{y^) = {q{y)E{y,0)yD{yD)q{0)) , 

S|(y2) = i {q{y)E{y, 0)yD\{0)) , (4.20) 
BUv^) = i (q{y)Eiy, 0)D{yD)q{0) ] 
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Expanding in y and averaging over y directions (Section |3.3|) we obtain the series 



Bl{y') = 



Bjy 
2D 



+ 0(/), 



Bliy'^ 



Bm 



D 



m 



2 2Bl + Bl + mB^ -m'^B^ 4 g 
){D + 2) 

/ + 0(2/'), 



B^ 2 Bl + 2Bl + 2mBl 



D 



12D{D + 2) 



BUy 



-y 



D" l2D{D + 2) 
We introduce 4 new bilocal functions 

Q'fiiy^) = i {q{y)E{y, 0)G^,(0)a^,g(0)) , 
Q'f2iy') = '-^^^ {q{y)E{y,0)Gx,,iO)axuqm , 

Q^h{y^) = ^ {q{y)E{y,Q)y^.G^Mluqm , 
Q'Uiy^) = ^ {q{y)E{y,0)y^G^,{OhMO)) . 

They are expressed via the basis bilocal condensates ( 4.2C1| ) as 
QVi(y') = -2B!{y') + 2m^B''fs{y'), 
Q'f2iy') = ^Bliy^) + mB^ {My') + 2y2/(,(y2 



2L> 



2D 



S|(y2)-4Dmi?3/^(y2)^ 



QV^iy') = ^ - ^l(y') ) - ^DmB^fsiy') - 2m'B-^fv{y') 



y^ 

Solving this linear system we obtain 



3 f„,2\ 



Bliy'^ 



Bliy') = ^ [Q'/2(y') - mQ3 (/^(y2) + 2y2/(,(y2 
= ^ [Q' (/4(y') - /3(y')) + 2m'Q^fy{y') 



bW^ 



y 

'2D 



Q*^Uy')+iDmQ^fs{y'' 



Substituting (|]2l|) and into ( p3| ) we arrive at the small expansions 

V + o(/), 



My') 


= 1 + 


f2{y') 


= 1 + 


h{y') 


= 1 + 


h(.y') 


= 1 - 



Ql-Ql-2Ql-m'Q^ 2 



2DQ^ 

2Ql-Ql-3Ql-mQ^ -m'^Q^ 2 



2(L' + 2)Q5 



6^ - 3Qf - 4Q| + 6Qi + 3mQl - Gm^Q^ 



+ 



12(D + 2)Q6 
eg? - 3Q8 _ QQl + 4Q^ 



12m (-Ql + Ql + 2Ql - mQ^ + m^Q^ 



24{D + 2)Qf 
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Now it is easy to expand the noncollinear quark condensate (|4.g| ) in small x keeping y fixed. 
We expand q{x), substitute x = {xy/y'^)y + x± (with x'j_ = A/y^), and perform the {D — 1)- 
dimensional averaging (Section |3.3| ) over x± directions orthogonal to y. Then we express the 
result in terms of the basis bilocal condensates ( [1.19 ), ( [4.22| ), and substitute ( 4.19 ), ( 4.24 ). 
Finally, we arrive at the result 



fs{x,y) = fsiy')-2ixy)f'siy^) + 



r/i(y' 



4(z)- i)g3y2- 

x'^y'^ - DjxyY / „, . 9. ^ 9 „„. o.\ m^A 
{D - l)y2 

fs{{y-xf) + 



2{D - l)y2 



/5(y2)+0(x^ 



A 



{D - l)y2 



m 



fT{x,y)=AD^{fs{y') + —fv{y 



2D' 



xy 

1/2 



Q 

■^$My') + f'siy') + V/^'(2/^) + ^ [My') + 2y^fUy 



m 



+ 0(x3), 
(4.26) 
+ 0(x'), 



g{y, x) = [h{y^) - Uy^)) + 6Q' [2{D + 2)My^) + Ay^ ^y^) + fv[y^)\ + 0(x), 



9{x,y) 



12Q3 



(Dfv{y^) + 2y2/y(y') " I^"i/5(y^ 



-^{g^(/4(y^)- 2/3(2/^ 
+ 2Q3 [6(1? + 2)/(.(/) + 122/2/v(y') - 4i?m/^(y2) + mVv(y')] } + 0(^2). 

If y is also small, we substitute the expansions of the bilocal condensates ( |4.2| ), ( [4.25| ) and 
obtain expansions consistent with (|4.9| ), ( 4.10| ), ( 4.12 ). This provides a strong check of our 
results. 

Expanding a noncollinear gluon condensate in x, we obtain terms like (Tr G^,y(y, 0)0^,^(0)) 
where G^,y{y, 0) = gG'^^E°-'^t^ and O^i/(0) is constructed from G^^(O) and D^. A straightforward 
generalization of the systematic classification procedure of Section to the case of such bilocal 
condensates seems impossible. Fortunately, it is not difficult to find all basis bilocal condensates 
up to dimension 6. Two dimension 4 condensates were introduced in Section [4.11: 



(4.27) 



A{y') = (T¥G^,(y,0)G^,(0)) , B{y') = (Tr GA^(y, O)Ga.(O)) 



There is one dimension 5 bilocal condensate 



C(y2) = (TrG^,(y,0)J^(0)) 



Dy, 



(4.28) 



and 5 dimension 6 ones 



Fi(y2) = i {TvG^,{y,0)G,x{0)Gxf.m , 



F2(y2) = i(TrG^,(y,0)G,p(0)G,^(0)) 



Dypya 



F^{y^) = i(TrG^p(y,0)G<,A(0)GA^(0)) 



Dypya 



i(TrGp,(y,0)G,A(0)GA.(0)) 



Dypya 



(4.29) 
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F4(y2) = (TrG^,(y,0)Z)^J,(0)) 
F5(y2) = (TrG^p(y,0)Z)^J,(0)) 



y2 



(the equivalence of the two definitions of F-i{y'^) follows from the C-conjugation symmetry). 
Small ip' expansions (^^), ([4 .61) of A(y^), B{y'^) were obtained in Section 4J; the other expan- 
sions are 



12(L> + 2) 



AD 



4(1) + 2) 
3G3_4 + 4G| — Gfi 



6 ,.2 



2(L> + 2) 



y' + o{y% 



F4iy^ 



F,{y^) = Gl + 



'-y' + 0{y'), 



2D 

W + o{y'). 



4(Ll + 2) 

Now we can derive the expansions in the bilocal gluon condensates 



DA 

y2 



Fx^^x^y ^yp 



(2Fi(y2) + F4(y2) + DA'{y^) + 2y^j^'{y^)) + 0(a 



2F3(y') + ^i^5(y') + + y^C'iy^) 



+ Z)i3'(2/2) + 2y^B"{y^) + 



_ B{y^ 



y' 



Fx^,x,y^x, = A G{y') + 2B'{y') - D 



My^)-B{y^) 



(4.30) 



+ 



y2 



SFsiy') - \F,{y^) - lc{y') - Sy'C'iy') 



^j^Aifl^^B^ + 2DA'iy^) - {3D + 2)B'{y^) - 6y^B"{y^) 



Fx^^M^^yu = A f 2S'(y2) _ DA'{y^) - D 
Axy 



+ 0(x^ 



A{y^) - B{y^ 

y2 



y' 



F,{y^) - 2F,{y^) - \F,{y^) - \c{y^) - y^C\y^) - 2Z?^i(^!L^^ 



+ 3DA'{y') - i3D + 2)B'{y') + 2Dy'A"{y') - 6y''B"{y'' 



Axy ( 2n n/iV 2x , . „// 2. ^r^ My"^) - Bjy'^) 
Fxf,,XuX^j.Xp = — 2- C{y ) - DA (y ) + AB [y ) - 2D ^ 



(4.31) 



{D + l)y4 



Z)A 3Fi(y2) + 2F4(y2 



+ i [xW -{D + 2){xyf) (-2F2(y2) + 6^3(^2) + F^{y^) + 3C{y^) + &y^C'{y^^ 
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+ ((£» + 4)x2y2 _ (5Z)2 ^ ^ 8)(.xy)2) 

- 2D (2xV -{D + 3){xyf) J^'iy") 



p,v ^puV ixXvV pX a 



A2 



2(i? + l)y2 



Solving the system of 6 linear equations, we finally obtain 



A 



Cl 



(D + l)(D-2)y2 



L>-3 



m(y') + F5(y') + C(y2) + 2y^C'{y') 

{D + 2)iDA'iy^) - 2B'{y^)) - 2y\DA"{y^) - 2B"{y')) + 0{x) 



1 



C2 



p + !)(/)- 2) 



(D - 5)(-DFi(y2) + 2F3(2/2)) - (D - l)F2{y' 



D-3 



DF^iy') + F^iy') + C(y^) + 2y2c"(y2) + d{D + 1) 



+ 2B'{y') - 2y\DA!'{y') - 2B"{y')) + 0(x) 



A{y^)-Biy^) 



C3 



xy 



+ 



D-2[y2 I 

1 x2 

2(L> + 1) y2 



D-3 



+ 



+ 2DF4y') - DF^{y') - (D + 2)(C(y^) + 2y'C\y')) 

- 2D(D + i\^kJ^_Eill + 2L>(D + 2)A'(y2) _ + 2L> + 2)B'(y'^) 

+ 4L>y2^"(y2) _ 4(D + 2)y^B"{y'') (4.32) 
1 (xy)2 



2L>(3Fi(y2) + 2F4(y2)) + + 2)(-2F2(y2) 

^(y2)-S(y2) 



2(Z) + 1) y4 

+ 6^3(2/') + i^5(/) + 3C(/) + 6y2c"(/)) + 8L>(D + 1) 
- 4D{3D + 4)A'(j/2) + 2(5^2 + lOD + 8)B'{y'^) 



4D{D + 3)y'A"{y') + 4(5Z^ + 8)y^B"{y'^' 



C4 - C5 



L>-2 



- C{y') - DA'iy') + p(F2(y') + Fs{y') 
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C4 + C5 



+ 



D-2 
1 xy 



y 

2D{D - 5)Fi(y2) + (1)2 _ 5)^^(y2) _ p2 ^^D- 22.)F^{y^ 



D + ly'^[ D-3 
+ 2Z)F4(y2) - 2F5(/) - p + 3)(C(y2) + 2y^C'iy^)) 



y 

.2 /|///„,2 



+ 2Z)(Z) + 3)y'A"{y') - 8{D + 2)y2i?"(y2) 



C6 



(I) + l)(D-2) 



(ZJ - 5){DF,{y^) -{D- mjy')) + {D - l)F2{y^ 

D-3 



+ DF^iy') -{D- l)F^{y') + D{A'{y') + 2y'A"{y')) + 0{x 



Substituting the y"^ expansions of the bilocal gluon condensates (|4.5|) , (4.6), ( 4.30 ) into ([4.3l| ), 
( 4.32| ), we obtain expansions consistent with ( 4.15| ), ( [4.16| ); this provides a strong check of our 
results. 



5 Heavy quarks 

5.1 Heavy quark condensates 

If the quark mass is large, all correlators can be expressed via gluon condensates only. This is 
also true for quark condensates (one-current correlators): 



This is an expansion in 1/m. 



Qk = ^Ckn{m)Gn. 



(5.1) 




Figure 4: One-loop diagram for a heavy-quark condensate 
In the one-loop approximation (Fig. ^) quark condensates are given by the formula 

' dp" 



Qk = {qOk[D^]q) = -i I [ (Tr Ofc[-ip^ - iA^]S{p)) 



(5.2) 



where Af^ is given by (O). We use the MS regularization: the space dimension is D = 4 — 2e, 



2^) D " V^tT ) (27r)^' 



(5.3) 
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/i is the normalization point, 7 is the Euler's constant. We substitute the quark propagator (|2.3|) , 
(|2.4| ) into (^]^) and average the integrand over p directions in the Z)-dimensional space. The 
result is expressed via the integrals 



i(-l)"(m2)2-" 



m 



2Ti)D{p^-m?Y (47r)2(n-l)! \ j ^^'^^"^ ^ + 



im / 1 



m 



(4^)2 \e 



--log^ + 1 +0(e), h 



Ir. 



+ 0(e) (n>2). 



(47r)2(n-l)(n-2)(m2)"-2 
Coefficients Cfc^ with d„ < contain ultraviolet divergences: 



I I ^ m? . I 

7fcn I - - log ) + 



(5.4) 



(5.5) 



We define the quark condensate QkifJ-) renormalized at the point n as the sum of the series (5.1) 
from which 1/e poles are removed. 

We obtain the following results |^ |6| 



1 



247r2 
1 



6m3(L + 1)N + —G^ ^(G? - 6G^) 

m 15m'^ 



+ ^(19Gf + 16G5 - 20G? - 78G| + 84G? - 4Gg - 18G?) + O ^ 

210m"' V?7i' 

^ 3mLG^ - —{Gf-Gl 



127r2 



m 



+ Y^i^Gt + 2Gl -Gl- + - - G?) + O ( — 



Ql 
Ql 
Ql 
Ql 
Ql 
Ql 

Ql 
Ql 



1 



127r2 
1 

~12^ 
1 

1 

4^ 
1 

~4^ 
1 



LGl + ^{l^Gl-\Gl-^G^,) + 

- 3mHL + l)G^ + -^Gf + O f 

2m / 

(G? + G8-4Gl) + of4 



mLGl + ^(4G8 - iG^ - G?) + O f 3 



(5.6) 



127r2 
1 



3m^(L + |)G^ + L{Gl - 2G^ - 2G^) + O ( ^ 



247r2 



12m2(L + 1)G? - L(G? - G^ - 2G^ + 2G| - 4Gi + G^) + O 



127r2 



Gg + O 



1 
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2 

where L = log As we have already mentioned, the series for the anomalous condensate 
A includes only d = 8 gluon condensates: 

A = -^{Gl + Gl-4Gl). (5.7) 

Many of these results can be obtained using various physical considerations § instead of 
the described "brute force" method. 



5.2 Heavy quark currents' correlators 



P 




Figure 5: One-loop diagram for a correlator 
Correlators of heavy quark currents can be expressed via gluon condensates: 



(5.8) 



The one-loop diagram (Fig. [S]) contains two propagators of the type (2^), (2^) (don't forget 
that one of them has the vacuum momenta sink on the other side, and is given by the formulae 
mirror symmetric to (^.3|), ( |2.4D !). After differentiations we obtain a formula with the integrals 

dk\ P{k) 

2^ In {k^ - m^Y'iik - p)^ 



2\m ■ 



£, [K- — fib- ]--\\K — vr — m 
The denominators can be combined using the Feynman's formula 

1 



(5.9) 



1 r(77. + m) 



n-1 



x'" "(1 — x)™ dx 
a"6™ r(n)r(m) J [xa + (1 - x)6]"+'" 



(5.10) 



After the shift of the integration momentum k ^ k + xp the denominator becomes k'^ + x(l — 
x)p'^ — m?; we may average the numerator P{k + xp) over k directions. The integrals reduce to 
the form ( |5.4| ) with — > — x{l — x)p'^ . We are left with one-dimensional integrals over the 
Feynman parameter x of the form Q 

1 

-- 



dx 



[1 +^x(l - x)f 



(5.11) 



(2n-3)!! 



n 



D! 



1 



2a 



^log 



V^ + 1 
V^-1 



n-l 



{k-l)\ 



2a 



n—k 



where ^ = —p'^/m?, a = 1 + 4/^. Using these formulae, one can calculate any gluon contribution 
to a heavy quark currents' correlator. The simplest example of d = 4 is considered in [|^; 
contributions with d = 6 and d = 8 were calculated in [01 . 
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6 Light quarks 



6.1 Limit m ^ in heavy quark correlators 

Let us consider the heavy-quark correlator (|5.8D (Fig. ^) at ^ m? . We can express it via 
both gluon and quark condensates: 

n(p) = nG(p) + nQ(p), (e.i) 

n k 

Here IIg{p) corresponds to the contribution of the region where the virtuahties of both quark 
and antiquark in Fig. ^ are large (fc^ ~ p"^), and I[q{p) corresponds to the contribution of 
the regions where either the quark or the antiquark has a small virtuality (fc^ ~ m?). These 
contributions are usually depicted as the diagrams of Fig. |[ 




Figure 6: Quark condensates' contribution to a correlator 
The quark condensates' contribution (Fig. ^ in the coordinate space is 

nQ(x) = {q{x)TS{x, O)rg(O)) + c.c, (6.2) 
where P is a 7 matrix, and c.c. means charge conjugate. In the momentum space 

^Qip) = (9(1 - i'D^do, - ^^a'Dpdadp + • • ■)TSip)Tq\ + c.c. (6.3) 



We substitute S{p) ( ^.3]) , ( |2.4D and average over p directions, and finally reduce (|6.3| ) to the 
basis quark condensates following the Section 3.1. The contribution of A is omitted because this 
condensate is a gluon one. The quark contributions with d = 7, d = 8 were obtained in 

Having obtained IIq{p), we can find Hcip) from the heavy-quark correlator n(p) using the 
expansions (|5.6|) : 

Ug{p) = ^ a'Jp^m)Gn = U{p)\g - nQ{p)\G 

n 

= an{p^,m)Gn - Y bk{p'^,m)ckn{'m)Gn. (6.4) 

n k,n 

The coefficients a„ (p^ , m) have singularities at m ^ arising from the regions k'^ ~ m? . These 
singularities have to cancel in (|6.4| ) giving a',^{p'^^m) finite at m ^ 0. Using the results of 
the gluon contribution with d = 8 was obtained in ||9|, |^ . 



6.2 Minimal subtraction of mass singularities 

The method of the previous Section is good when the heavy quark correlator is already known. 
But when we want to calculate a light quark correlator from scratch, there should be a simper 
way than to solve a more difficult heavy quark problem first. Such a method was proposed in 
and used in for d = 8 calculations. 
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Now we want to go to the limit m — > before D ^ 4. The difference of a renormahzed 
quark condensate and a bare one is (see (|5.5|)) 



Qk-Qt'''' = -^ E "i'''"''"7fcnG„ ^ Yl ^kuGn. (6.5) 

dn<dh d„=dk 

At m = all loop integrals for Q]^'^'^^ vanish because they contain no scale (ultraviolet and 
infrared divergences cancel each other). Therefore we obtain for the gluon contribution to a 
correlator (see (|6.4p) 

UG{p)=Uip)\G + ^ J2 hiphknGn- (6.6) 

dn=dk 

Here n(p)|G' is calculated in the MS scheme with m = 0, and 7fc„ are mixing coefficients of quark 
condensates Qk with gluon condensates G„ of the same dimension (coefficients at L in (5.6)). 
Omitting 1/e poles we finally obtain 

UG{p)=Uip)\G+ Yl (6.7) 

dn=dk 

The quark condensates' coefficient functions 6^ should be calculated at m = up to linear 
terms in e (using D-dimensional tensor and 7 matrix algebra, and in particular D-dimensional 
averaging, see Section ^ ). 
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